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Abstract 

Lattice artifacts in the 2d 0(n) non- linear cr-model are expected to be of the form 
0(a 2 ), and hence it was (when first observed) disturbing that some quantities in 
the 0(3) model with various actions show parametrically stronger cutoff depen- 
dence, apparently 0(a), up to very large correlation lengths. In a previous letter [1] 
we described the solution to this puzzle. Based on the conventional framework of 
Symanzik's effective action, we showed that there are logarithmic corrections to the 
0(a 2 ) artifacts which are especially large (In 3 a) for n = 3 and that such artifacts 
are consistent with the data. In this paper we supply the technical details of this 
computation. Results of Monte Carlo simulations using various lattice actions for 
0(3) and 0(4) are also presented. 



1 Introduction 



In a previous letter pQ we presented results on logarithmic corrections to 0(a 2 ) 
lattice artifacts for a class of lattice actions for the non- linear O(n) sigma- model in 
two dimensions. It is the purpose of this paper to supply the technical details of 
the computation. The main results, which are summarized in subsection 3.5, are 
that the generic leading artifacts are of the form a? [ln(a 2 )]"^ n , and this result 
together with the next-to-leading expressions describe well the lattice artifacts in the 
step scaling function [1] , which are for n = 3 in a large range of the cutoff apparently 
of the form 0(a) [2j. The goal of our work was indeed to explain this long-standing 
scientific puzzle which was mentioned by Hasenfratz in his lattice plenary talk in 
2001 [3]. 

Most of our knowledge concerning renormalization of quantum field theories 
stems from perturbation theory. Although there are no rigorous proofs in gen- 
eral, many of the results are structural and hence considered to carry over to non- 
perturbative formulations. Indeed there is supporting evidence from various studies, 
e.g. of soluble models in 2 dimensions and of 1/n expansions of some theories. The 
same situation holds concerning cutoff artifacts in lattice regularized theories. 

Thus we start our paper with perturbative considerations. In section 2 we 
briefly summarize perturbative renormalization of the sigma model in the framework 
of dimensional regularization, including the renormalization of isoscalar composite 
operators of dimension 4. Next we discuss a large class of lattice regularizations. 
We summarize the results for connected 2- and 4-point functions to 1-loop order 
and also the relations of the corresponding renormalized functions to those of the 
dimensionally regularized theory. 

In the early 80's Symanzik was working on the nature of lattice artifacts, in 
particular with respect to his improvement program |4|5|6|7|8] . In this paper the 
general theory is not discussed; we only consider in section 3 Symanzik's theory 
applied to the 2-dimensional O(n) cr-mode0. Nevertheless the spirit of the gen- 
eral theory can already be understood by studying this example. Symanzik's main 
conclusion is that leading artifacts are summarized in an effective action. 

In this framework generic lattice artifacts are, in particular for asymptoti- 
cally free (or trivial) theories, expected to be integer powers in the lattice spacing 
0(aP) ,p = 1, 2, . . . up to possible multiplicative logarithmic corrections. In partic- 
ular this framework explains why in the 2-dimensional cr-model quadratic artifacts 
are expected. To obtain the effective action in lowest orders of perturbation theory 
we compute the leading lattice artifacts in the 2- and 4-point functions to 1-loop 
order and express them in terms of insertions of composite operators of dimension 

Concerning Symanzik's program for QCD, see e.g. [9] and references therein. 
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4 in the continuum. This is a rather lengthy computation and thus many technical 
details are relegated to the Appendices. 

In section 4 we reanalyze the presently available Monte Carlo data on the step 
scaling function for cases n = 3 and n = 4. In particular for the case n = 3 we study 
three different lattice actions, and see how the particular logarithmic corrections 
predicted by Symanzik's theory solves the puzzle of apparent linear artifacts. Similar 
conclusions that logarithmic corrections may explain the puzzle have been reached 
previously in studies of the O(n) model in the first orders of the 1/n expansion 
|10|11|12] . In Appendix E we describe a modified version of Hasenbusch's improved 
estimator which was essential to obtain a sufficiently small error for small a/L values. 

There is a vast literature on the subject of this paper, however we feel that 
this technical report is not the appropriate place to properly review all important 
contributions. Hence we restrict ourselves here to citing a selection of papers where 
the reader can find further references. 2-dimensional O(n) spin systems, including 
the Ising model and the large n limit, are reviewed in [13] . Various aspects (on the 
lattice and in the continuum) of the exactly solvable n — > oo limit, using a number of 
different lattice actions, are reviewed in [14]. Symanzik's improvement program for 
the 2-dimensional O(ra) models are studied in |15|16j and some exact results about 
the lattice artifacts in the large n limit are given in [12J. 

2 The non- linear 0(n) sigma model in two dimensions 

The non-linear O(n) sigma model describes the interaction of spin fields S a (x) , a = 
l,...,n with the constraint S(x) 2 = 1 in two space-time dimensions. Here we 
consider the space-time to be Euclidean. Formally the Lagrangian is that of a 
free massless field and the interaction derives entirely from the constraint. The field 
theory requires an ultra-violet regularization; in the next subsection we will consider 
dimensional regularization and then in subsect. 2.2 a class of lattice regularizations. 

2.1 Dimensional regularization 

The action of the O(n) model in D = 2 — e dimensions (including source terms) is 

A=^Jd D x!^d l ,S-d l ,S-I-Sj, (2.1) 

where go is the bare coupling. To satisfy the constraint it is usual to parameterize 
the bare spin field by 

S*=Sb7r\ i = l,2,...,n-l; S n = * = y/l- g 2 ^ . (2.2) 
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The source dependent action (|2.ip appears in the generating functional 

Z[I}= [ (Vn)e~ A , (2.3) 



which can be used to obtain bare correlation functions of the field S a (x) using the 
formula 

""*'"' ^l-^Wjj^j- w4jL^'' (2 ' 4) 

where the functional derivative is taken at 

I$(x) = m 2 5 an , (2.5) 

i.e. a mass term (external magnetic field) is introduced to avoid infrared singularities. 
As proven by David [T7], O(n) invariant correlators are infrared finite order by order 
in perturbation theory and hence for these the limit m — > can be taken at the end 
of the calculation. 

In their seminal paper [18] Brezin, Zinn-Justin and Le Guillou prove the renor- 
malizability of the O(n) model using functional methods. In this paper they show 
that the generating functional iJ -1 [/o]>Z[^] is finite as function of the renormalized 
quantities j a (x),g, fj,, tur if we write 

I a (x) = Z 1 (g,e)Z- 1 / 2 (g,e)g 2 f(x)+lS, 

g 2 = ^Z 1 (g,e)g 2 , (2.6) 
m 2 = Z 1 (g,e)Z-^ 2 (g,e)m 2 R , 
where the renormalization constants contain only pole terms, 



Zi(g, e) = 1 + — 2~ + • • • > 

see 1 

7( v _ , 705 2 7i9 4 , g 4 (7o + 2/3q7o) , 



(2.7) 



with 



A = V i . ft = ^ (n + 2Kn-2) 

4 \ Pl 8vr2 ' ^ 2 64tt3 ' (2.8) 

n- 1 v ; 

7o = — , 7i = , 

where the 3-loop coefficient /?2 will appear in (|2.13p below. Functional derivation 
with respect to the source j a {x) gives renormalized correlation functions, i.e. cor- 
relation functions of the renormalized fields S R = Z~ x l 2 S a . The relation between 
bare and renormalized correlation functions is given by 

G?B) (9, A*) e) = Z- r/2 (g, e)G x (g ,e) , (2.9) 
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where the upper index X symbolizes any r-point correlation function (in x-space or 
in Fourier space). We assume that X is O(n) invariant and that the m — > limit 
has been taken. Finiteness means that the limit 

Gf R) (g,fi) = gf R) (g,(i,0) (2.10) 

exists and defines the renormalized correlation function in two dimensions. 

The renormalization group (RG) equations express the fact that the bare cor- 
relation functions are independent of the renormalization scale /i. In terms of the 
renormalized correlation functions this is expressed as 

{D + ^(g)}gf R) (g,»)=0, (2.11) 

where the RG differential operator is 

2, ="| +/J(9) I' (2 ' 12) 

and the RG beta and gamma functions are defined as 

m = f - - £ d U (9 , e) = -/w - ^ -P2 9 7 +... (2.13) 

Z + 9 dg 

and 

7(5) = {m - f } = W + W + ... (2.14) 

We also introduce the RG invariant A-parameter in the MS scheme by the formula 



where 



and 



A M - s = eiA MS , A MS = t* {WoQ 2 ) X e ^ S 9 \ (2.15) 
X = ^ = ^~ r 7 = ln4 ?T + r(l), (2.16) 

fMf W ^} (2 ' 17) 
In the O(n) model, instead of Aj^g we could also use M, the physical mass of the 
particles, since its relation to the Lambda parameter is known |19| : 

M /8\ x 1 

(2.18) 



a ms W r(i + x) 
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2.1.1 Local operators 

We now turn to the renormalization of local operators. In our Euclidean framework 
operators Oi correspond to insertions into the generating functional: 

A 



O t — > Zi[I]= I (pir) Oie-* (2.19) 

Beyond operators corresponding to local expressions of the O(n) field S a (x) and its 
derivatives we will also consider here operators Oi depending on the source I a (x). 

Bare correlation functions with Oi operator insertion can be obtained by using 
Zi[I] in (|2.4p in place of Z[I\. We will denote these symbolically as . 

Operators are mixed with other operators under renormalization. Renormalized 
operators O^m-s are given by 

O i{R) =Z ij (g,e)O j , (2.20) 

which is a symbolical expression of the fact that 

Z- 1 [I ]Z ij (g,e)Z J [I] (2.21) 

is finite in terms of j a (x),g,fi and m R . In terms of correlation functions with 
operator insertion we have 

G? {R) (g,»,e) = Z- r / 2 (g,e)Z ij (g,e)gf(g ,e) (2.22) 

and finiteness means that the limiting correlation functions 

G?(R)(9,V<) = Gf {R) (9,»,Q) (2-23) 
exist. They satisfy the RG equation 

{v + ^(g)}gX {R) (g,ri + ViMG?(R)(9,») = 0, (2.24) 

where the anomalous dimension matrix is defined by 

eg\ dW sj (g,e) 



Ui j (g) = Z is (g,e)[(3(g)-f) 1^ (2.25) 



with Wij the matrix inverse of Zij. 



In perturbation theory (PT) we have 

g 2 g 4 g 4 

Zij(g, e) = Sij - —Wij + —Pij + ^ ( w isWsj + ZPoWij) + • • • (2.26) 



and 

nj{g) = -Wij9 2 + Pij9 4 + • • • (2.27) 
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2.1.2 Dimension 4 operators 

Brezin et al. prove in [18] that the following set of mass dimension four, O(n) 
invariant, Lorentz scalar operators is closed under renormalization. 

1 = ~(d lJl s-d tl s) 2 , 

2 = \ (d p s ■ d u s) (d p s ■ d u s) , 

o 

3 = ^as-as, (2.28) 

O^-ad^S-O^S, 

where 

a = 9z±HMl. (2 . 29) 

a 

Although the source dependent terms O4 and O5 look O(n) non- invariant it is 
demonstrated in Appendix A that in fact all (otherwise O(n) invariant) correla- 
tors containing an insertion of these operators are O(n) invariant. Brezin et al. 
[18j showed that they must be included in the operator renormalization scheme for 
consistency for m/0 and off shell. The one-loop mixing matrix for dimension four 
invariant Lorentz scalars is calculated in the paper using this basis. This result and 
further considerations on the 5x5 operator renormalization problem are discussed 
in Appendix A. 

In our work we also need to consider four-index symmetric tensor operators: 



t^p, = S • dpdvdpdaS , (2.30) 

k„u P a = l{(d IJt s-d u s) (d p s-d a s) 

+ (d p s ■ d p s) {d v s ■ d a s) + (d p s ■ d a s) (d v s ■ d p s) } . (2.31) 

In particular in the Symanzik effective Lagrangian we will encounter the following 
operators A, B which are defined from the totally traceless parts t, k of t, k: 

D 

A = ' (2.32) 

D 

B = Y$m*' ( 2 - 33 ) 
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Note that these operators are invariant only under discrete (D-dimensional lattice) 
rotations. Their precise definitions and a discussion of their renormalization are 
given in subsection (A. 2). 



2.2 Lattice regularization 

In the lattice regularization the fields are restricted to the sites of a regular square 
lattice. With the usual assumption of universality an infinite class of local lattice 
actions could be invoked. In this paper we will only consider O(n) lattice actions 
quadratic in the spins: 

A = f fl4 E E S\x)K(x - y)S\y) , (2.34) 
x,y b 

with K short range and demanding 

Y,K(x) = 0, (2.35) 

X 

and 

K{z) = K(Kz) , (2.36) 

where 1Z is a lattice rotation or reflection. 

Let K p be the Fourier transform of K(x): 

K p = a 2 Y,e~ ipXK ( x )- (2-37) 



X 



The behavior of K p for small a is assumed to take the form 

K p = p 2 [1 + a 2 r(p) + 0(a 4 )] , (2.38) 

with 

r(p) = p 2 [K!R(p) + k 2 ] , (2.39) 

4 

R{p) = , (2.40) 

where we have introduced the notation p r = X^P^- 
For the familiar case of the standard action (ST) 

K(z) = K ST (z) = ^ [2**,o " S Ztafi - Sz-a/i] , (2.41) 
ft 

where /t is the unit vector in the \x direction, with Fourier transform 

K p; ST=p\ ^ = ^sin(^), (2.42) 
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and so in this case k\ = — 1/12 , «2 = 0. 

Expectation values of an arbitrary Euclidean observable 0[S] are given by 



<0> = i J[dS]e- A 0[S], (2.43) 



where 



[d5]=IJdM(5(x)). (2.44) 

Here dfi(S) denotes the O(n) invariant single spin distribution normalized to 1, 

d/x(S) = 5 (S 2 - 1) [J dS a (2.45) 

a 

and the partition function Z is such that (1) = 1. 

Perturbation theory is an expansion for (3 — ► oo and so for this purpose we set 

0=^2- (2-46) 

As in the case of dimensional regularization, to avoid intermediate divergent ex- 
pressions (from some Feynman diagram contributions) it is useful to introduce an 
infrared regulator e.g. work in finite volume or add a coupling of the spins to an 
external magnetic field. These aspects are discussed further in Appendix B. In this 
section we give results for infinite volume and zero magnetic field. 

In our work we consider correlation functions of scalar products of the spin field, 
and for this purpose it is convenient to define 

9(x,y) = S(x)-S(y)-l. (2.47) 
2.2.1 One 6 correlation functions 

One 6 correlation functions have a perturbative expansion of the form: 

C(x,y) = {6(x,y)) (2.48) 
= (n-l)X 2 J2^ r Cr(x,y). (2.49) 

r=0 

In the lowest order (for p ^ 0) the Fourier transform is given by 

C (p) = ^ = ^[l-a 2 r(p) + ...] . (2.50) 
j\ p p 



In the next order we have 

C l (p) = K- 1 



(2.51) 
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where F(p) corresponds to an eye diagram 

F(p) 

and k p to a tadpole 



K p — K s — K p+S 
Kp+sKs 



k p = J ^(K p + K s -K, 



p+s) 1 



where we have introduced the shorthand notation 

w/a ^2 g 



it /a 



(2vr) 2 



(2.52) 
(2.53) 
(2.54) 



(2.55) 
(2.56) 

(2.57) 



2.2.2 Two 6 correlation functions 
Consider the ^-connected correlation function 

C{x 1 ,yi;x 2 ,y 2 ) = (0(x 1 ,y 1 )6{x 2 ,y 2 )} - (9(x 1 ,y 1 ))(e(x 2 ,y 2 )) ■ 
This has a perturbative expansion of the form: 

C(x 1 ,y 1 ;x 2 ,y2) = (n - 1)Aq ^ \l r C r {x u y x \ x 2 , y 2 ) ■ 

r=0 

Define its Fourier transform by 

C r (pi,qi; P 2,q2)=a 6 £ e -^ + ^ + ^ + ^C r ( Xl , Vl ; x 2 , y 2 ) . 

To avoid disconnected contributions we restrict attention to momentum configura- 
tions such that 

Pi + qi + P2 + <?2 = , 

pi^O, P2 ^0, qi ^0, q 2 ^0, 

p 1+P2 ^0, Pl+q2 ^0, (2.58) 

i.e. all momenta unequal zero and the sum of pairs of momenta associated with 
different 0's unequal to zero. Then 

C (pi,qi;p 2 ,q 2 )=0. (2.59) 

In the next order we have simply 



Ci(pi,qr,p 2 ,q 2 ) 



1 



1 



K K K K 



K, 



1 1 
+ 



P1+P2 



K, 



Pl +92 



K K K K 



[(n - l)K pi+qi + K P1+P2 + K Pl+q2 ] . (2.60) 



The expression for the next order which is rather lengthy is given in Appendix B. 
The continuum limits of these results and the leading lattice artifacts will also be 
considered there. 
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2.3 Relation between the renormalization schemes 



Renormalized lattice correlation functions which have a finite continuum limit (a — ► 
0) can be obtained by applying wave function and coupling renormalization. In 
particular the correlation functions in the MS scheme are obtained by 



C m (pi,qr, ■ ■ ■ ;Pk,1k) = Z n K C{p 1 ,q 1 \ . . .;p k ,q k ) + 0(a 2 



(2.61) 



provided one expresses the lattice functions in terms of the renormalized MS coupling 
related to the bare lattice coupling through 



„2 _ V-U2 
9 — ^i R A o ■ 



(2.62) 



We of course checked explicitly that our computations confirmed this claim. 
The renormalization constants have a perturbative expansion of the form 



Z R 



i + z R i y + 4V + ..., 

i+z&y+*£V+.... 



At 1-loop one obtains: 

^R 



z 



(1) 

1R 



(n - 1) — 
(n-2) — 



ln(a/i) + -(ci +7) 
ln(a/x) + - (ci +7) 



(2.63) 
(2.64) 

(2.65) 
(2.66) 



where the constants Aq,c\ are given in (irTT7|) . (lrl22l respectively. Defining the 
lattice A-parameter by 



A latt = a- 1 (2^\l)~^ [1 + 0(A§)] , 



one gets the relation 



A 



latt 



Ams 



exp 



i(c 1 + 7 )-^ 

2 y " n-2 



(2.67) 



(2.68) 



For the standard action this result was first obtained by Parisi 



3 Symanzik's theory of lattice artifacts 
3.1 The effective Lagrangian 

We write the lattice Lagrangian including the source terms symbolically as 



-latt — -^y2 ' C^) 1 ^ 1 



JS. 



S 2 



1. 



(3.1) 
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with some lattice regularization of the kinetic term. We use the corresponding action 
in the generating functional for bare lattice S'-field correlation functions, which, after 
Fourier transformation, become functions of the momenta, the bare lattice coupling 
Ao and the lattice spacing a. Performing, for fixed momenta and coupling, a small 
a expansion we can write 

QL(M = S X(0) (A ,a) +a 2 S x «(A ,a) + (a 4 ) , (3.2) 

where both the scaling functions Q x (°> and the leading cutoff corrections Q X W are 
still weakly (logarithmically) depending on a. 

The separation of the full lattice correlation function into a scaling piece and 
cutoff corrections is unambiguous and straightforward in PT. In fact, in PT at I 
loop order both terms are finite (order t) polynomials in In a. For example (using 
the results in Appendix B) for the 2-spin correlation function at 1-loop one has 

Ci{p) = 7—5 \L{ap) + ci + 4vrA 
4-717^ L 

+ aV Tt - R (p)j L M + C 3 R (P) + C4 j + O(aV)] , (3.3) 
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where 



L(ap) = ln(aV) , (3-4) 



and all undefined constants are given in Appendix B. 

One of our main assumptions here is that the expansion (|3.2p makes sense also 
beyond PT. Usual renormalization theory deals with the scaling part Q x \°) and in 
later subsections we will analyze the next term Q x v-> using Symanzik's method. 

Symanzik's local effective Lagrangian C e a which is defined in the continuum in 
D dimensions, and is written in terms of the bare S'-field normalized to unity, de- 
pends nonlinearly on the source I and has the property that the generating functional 
obtained from it completely reproduces the correlation functions corresponding to 
(I3.ip up to terms of order a 4 . It is of the form 

— £cS = ~7T~o {dnS ■ daS) H tjI ■ S 

,2 



+ ^ [Y l 1 + Y 2 2 + Y 3 3 + Y A A + Y B B} 



3 - - ^ 

+ 

9o 



-2 { w x i ■ us + w 2 i ■ s (d^s • d^s) } 



a 

7io 



+ - { XJ 2 + 



X 2 {I-Sf) 
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Here the operators 0\, O2 and O3 are defined in (|2.28j) . the operators A and B 
in (|2.32p . (|2,33p and the constants Y\,...,X.2 can (in principle) be determined by 
comparing correlation functions calculated directly on the lattice and from (j3.5f) . 

Actually, this form of the local effective action for the O(n) nonlinear sigma 
model cannot be found in Symanzik's published papers |4|5j . But in the case of the 
4 model an analogous local effective Lagrangian is introduced in [1] and in [5] an 
improved lattice action is presented for the O(n) sigma model. (13.5j) is constructed 
using the cp 4 model analogy and has exactly the same terms as the improved action 
in [5]. It is not unique: as discussed in [5] there is an ambiguity in choosing the 
coefficients which can be used to put (for example) W2 and X2 equal to zero. To 
be able to identify the lattice correlation functions with the ones obtained using the 
effective action we also have to relate the couplings Ao and g and the sources J a and 
j a , as will be discussed below. We can eliminate the terms quadratic in the source 
fields by using the identities derived in Appendix A. We get 



1 1 2 f 5 

£cff = ^ S -d^S) + -I -S + ^l Y A A + Y B B + YYiOi 

l % % So { i=1 



(3-6) 



where 



Y A = Y A , 
Y B =Y B , 



Y x = Yl + 8W 2 + 8X 2 , 
Y 2 = Y 2 , 

F 3 = % + 2Xi - 2Wi , (3.7) 
F 4 = 2W 2 + 4X 2 + 4Xi - 2Wi , 
F 5 = 8X 2 + 8X1 . 
One can see that Y \ are invariant under the transformations 

SX X = -w, 5X 2 = to, 5Y X = -8to , 5Y 3 = 2uj (3.8) 

and 

SWi = p, 5W 2 = p, SY\ = -8p , 5Y 3 = 2p . (3.9) 

These correspond to the ambiguities of the original effective Lagrangian (|3.5p dis- 
cussed above. 

Finally we introduce the operator basis consisting of 

Ui = ^Ui, i = l,...,5, (3.10) 
% 

together with Uq and U7, where these and the operators C7j ,i = 1, . . . , 5, are defined 
in Appendix A. In this basis we have 

1 1 7 

- £ eff = --^ (d^S ■ d^S) + -I.S + a 2 Y, Yfa, e)Ui , (3.11) 
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where the coefficients Yi are appropriate linear combinations of Y i and Y a,b- This 
particular basis of operators is chosen such that they are renormalized according to 

U t{R) =K ij {g,e)U 3 , (3.12) 

where the matrix of renormalization constants is block diagonal consisting of the 
5x5 block for i,j = 1, 2, 3, 4, 5 discussed in subsection (A.l) and the 2x2 block for 
i,j = 6,7 discussed in subsection (A. 2). Recalling the redefinition (|3. 10|) we have 
for the first block 

K ij = Z 1 K ij , (3.13) 

where K is the renormalization matrix for the U{. 
We now define 

7 

c j (g,e)=J2Yi(9,£)K^(g,e) (3.14) 
i=i 

and using this definition we can write 

7 7 



i=l i=l 

which shows that the limit 



must exist. 



J2Y i U i = ^c i U i(R) , (3.15) 
i=l 

c i {g) = c i (<7,0) (3.16) 



3.2 Relations between correlation functions 

We now write down the precise relation between the lattice correlation functions and 
the ones obtained by using the effective action. We identify (for simplicity) the scale 
parameter /i of dimensional regularization with the inverse of the lattice spacing a. 
In this case we can write 

gX tt (\ ,a)=y r (g)gf R) (g,a- 1 ) 

7 (3 17) 

i=l 

where the finite wave function renormalization constant y(g) comes from the rela- 
tion between the lattice source J and the (renormalized) dimensional regularization 
source j: 

f{x)=y{g)J a (x) (3.18) 

and from universality it follows that a relation of the form g = G{\q) must exist 
between the couplings of the theory. 
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Comparing (|3.17p with the expansion (|3.2p we see that 

g X{0) (X ,a)=y r ( g )gf R) (g,a- 1 ), 
7 



(3.19) 



i=l 

This result can also be written as 



Q&*(*o,a) = G X{0) (Xo,a) {1 + a 2 8 x (X ,a)} + O (a 4 ) , (3.20) 

where 

7 

S x (X ,a)=J2c l (g)Sf(g,a) (3.21) 
i=i 

and 

3.3 Tree level effective action 

Let us write down the tree level effective action. If we start from the standard lattice 
regularization of the sigma model with action 

- a st = i £ (s. ■ s I+ , - 1) — (^V^) 2 ■ < 3 ' 23 > 

we get classically in the continuum 

" -4st = -i £ /d 2 x (d,S) 2 + -^^/d 2 x(S. a 4 5) + O (a 4 ) . (3.24) 

Written in the basis of operators introduced in subsection (2.1) this corresponds to 

- & = < 8 ^ • *" s > + 1 {h A + + M • < 325 > 

If we start from a different lattice action in our class (which is still quadratic in 
S), we will find in the tree level effective Lagrangian a different linear combination 
of the two operators A and O3, including the possibility of the vanishing of both 
coefficients (which occurs for example for the improved action). In general we have 



42 = • VO + I {£» + >} + O M , (3.26) 
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where 

e = e 4 = 1 (3.27) 
for the standard action ST. For the general case 



e = -12KJ , e 4 = — 12/si - 16k 2 • (3.28) 

tie leading coe 
is the leading term in the expansion 



We can read off the leading coefficients c- ^ of the operators in our basis. This 



1=0 

Using (the third row of) (|A.30|) and (|3.26[) we find 

(o) _ e 
Cl ~ 4(n - 1) ' 
(o) _ e 



(5)=E4V- (3.29) 



.(0) = e£ 
24 

0. 



(3.30) 



2 4(n - 1) ' 

(0) _ e (n 2 - 2n - 4) 

c 3 - 4n(n-2) ' ^q) 

(0) 9ep 7 

° 4 ~ 4(2n - 1) ' 

(o) _ (3 - n)(3 + n)e 

Cs ~~ 16n(n - 1) 

Note that cf ] = for n = 3. 

3.4 Renormalization group considerations 

In this section we will analyze the structure of the lattice artifacts using RG methods. 
We recall that operators are renormalized in our basis according to (|3.12[) where the 
operator renormalization matrix is of the form 

2 4 4 

Kij(g, e) = Sij - — kij + — vf^ + (k is k s j + 2(3 kij) + ... (3.31) 

and the anomalous dimension matrix is defined as 

*M = K is(g, e) (m - f ) = -h 3 g* + v^g* + vfg« + ..., (3.32) 

where My is the matrix inverse of Kij. It is clear that if the matrix K is (block) 
diagonal, then so is the matrix v and if the matrix K is (block) triangular, then so 
is the matrix v. 
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kij (in (|3.32p ) is a diagonal matrix of the form 

% = 2p \S tj (3.33) 

with eigenvalues corresponding to 

I ft 1 — Tl I 1 

A,:= ^;-l;0; -; -;0;-l . (3.34) 

[n — 2 n — 2 n — 2 J 

It is easy to see that the functions 5f (defined by (|3,22p ) satisfy the RG equation 

VSf (g, a) = -Vijig) Sf(g, a) , (3.35) 

where 

To solve this partial differential equation we introduce the matrix Uij(g), which 
solves the ordinary differential equation 

UUg) = -p ls (g)U sj (g), (3.37) 



where 

This has the expansion 
where 



2Aj 

Pijia) = — - Sij + Rij{g) , (3.39) 



9 



oo 



%(s)=£/f S^ 3 - (3-40) 



«=2 



The two-loop and three-loop coefficients are: 

1 (2)_/?l 



{2) --— u^-^k- 
Pij ~ a V ij R 2 ^3 ' 



(3.41) 



9x3 Po v + /3 2 * U $ 
If we find the solution of (|3.37|) we can write the general solution of (|3.35|) as 

5f(g,a)=U ij (g)Df(A), (3.42) 

and the lattice artifacts are of the form 

7 

5 x (\ ,a) = Y J M9)D?(A), (3.43) 



8=1 
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where 

7 

v i (g)=J2c s (g)U si (g). (3.44) 

s=l 

The functions Dj depend only on A, the RG invariant combination of g and a. These 
functions are non-perturbative and depend on the quantity (X) we are considering. 
On the other hand, the coefficients ii(g) are perturbative and universal in the sense 
that they remain the same for all physical quantities (but depend on the lattice 
action we started with). 

We take the following ansatz: 

Uijig) = {<% + Qijig)} g- 2A > , (3.45) 
where Qij(g) has the loop expansion 

oo 

QM = T, k< S ) 9 2l ~ 2 . (3.46) 

1=2 

Here the coefficients kf^ may still weakly (logarithmically) depend on the coupling. 
This can arise if the difference between two eigenvalues Aj — Aj is a non-zero integer, 
which is possible in our case, i.e. for n = 3. We can however ignore this subtlety, 
because we verify in Appendix D that for the quantities we need here it plays no 
role. 

We can now write the lattice artifacts as 

7 

5 x (X ,a) = Y J v l {g)g- 2 ^ D?(A) , (3.47) 
i=i 

where 

oo 

Vi(g) = ci(g) + c.(g)QM = £ vf ] g 21 . (3.48) 

The spectrum of one-loop eigenvalues given by (|3.34p plays a crucial role in our 
considerations. The leading term corresponds to 

Tl 

Ai = = n X = 1 + 2 X , (3.49) 

n — 2 

and the sub-leading one to 

A 5 = ^— = X - (3.50) 
n — 2 

We thus have the leading expansion 

5 x (\ ,a)= Vl D? {g- 2 ) l+2x + v,D x {g- 2 f + ... 

= D? {(g- 2 ) 1+2x v^ + vP(g- 2 ) 2 x} + 0((g-r) . ^ 
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We have already noted that for the n=3 case = = 0. This means that the 
corrections start one power later and in this case we have the leading expansion 

5 X (\ , a) = D* {vf g-« + v? g~* + <T 2 } + 0(1) . (3.52) 

The expansion coefficients are 

„(o) _ jo) _ eo 

Vl ~° l ~4(n-l)' 

v i - c i + c * K si i (3.53) 



Concretely we have 



= 2(A,-l-A.) ^ • P - 54) 

which is different from zero for s = 1, 2 only. 

We first note that the functional form (|3.51|) . (|3.52p is only valid for actions 
where the leading coefficient = q does not vanish. An important special case 
is Symanzik's (tree) improved action, where the above condition is not satisfied. We 
also note that, as can be seen from (|A.12jl . the connected correlation functions of 
the operator U§ are pure contact terms and therefore do not contribute to on-shell 
physical quantities. For such observables the sub-subleading corrections in (|3.51|) 
come from the operators U3 and Uq and are 0(1) (for n > 4). 

We also need the connection between the lattice coupling Ao and g: 

g 2 = Xl + ^Xt + ..., (3.55) 
which can be calculated from (|2.66p by setting a/x = 1. We find 

if] = --( Cl +j)(n-2) + 2irA . (3.56) 
We will also use the inverse coupling = 2-k/Xq. 
3.5 The final result 

The information above is all we need to write down the final result: 



* X < V, ) = _ Z?(A) j (p) WX + rW (ff t \ + O ( ^ ) (3.57, 
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for n > 40, and 

S X (X , a) = ^3 D±(A) + r^p 2 + r< 3 >£} + 0(1) (3.58) 
for n = 3. We write the two-loop coefficient as a sum of three terms: 

r (2) = r (2) + r (2) + r (2) ^ (3 ^ 5g) 



with 



- 8 " (n " 1} -c« , (3.60) 



.(2) 



11 n -2 

.(2) _ ,„ x 2 



eo 

(1-V), (3.61) 



(2tt)< 



n — l n 



(3.62) 



The above form (13.57p . (|3.58p of artifacts is general in the sense that it holds for all 
observables: the leading perturbative coefficients are the same for all observables 
(but depend on the lattice action chosen) . Only the overall coefficient of this univer- 
sal series depends on the quantity in question. This coefficient is nonperturbative 
and for any physical quantity X it can be parameterized by 

D?(A) = A 2 f x (z), (3.63) 

where z = LA or a vector with components zj = LjA, where L (or the Lj-s) are the 
relevant length scale(s) in the problem. 

If we use the RG relation between the lattice spacing a and the inverse coupling, 
we can write the functional form of the leading artifacts as 

a 2 [Ha 2 )] n/{n - 2) . (3.64) 

Note that our final result is completely consistent with the results of refs. |10|11|12] . 
who found leading artifacts oc a 2 ln(a 2 ) in the leading and next-to-leading orders of 
the 1/n expansion. 

We have calculated the three terms contributing to the sub-leading coefficient 
(|3.59p . For the calculation of the first term we need the 1-loop coefficients of the 
effective action. We obtained these using the computations described in section 2 
of the 2- and 4-point spin correlation functions both on the lattice and in the con- 
tinuum. More precisely we need the scaling part and the 0(a 2 ) piece of the lattice 
correlation functions and in the continuum we need the original correlation func- 
tions as well as the ones where those dimension four operators that appear in the 



For on-shell physical quantities the corrections are O(l) only. 
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tree level effective action are inserted. These are both long computations and are 
discussed in more detail in Appendices B and C respectively. 
For the spin-four tensor operators we find 



c 6 - 2 Al ~ ~2 7 ~ 16 17 ' 



The definition and numerical values of these and all other lattice integrals can be 
found in Appendix B. 

The coefficients of the scalar operators are simplest in the original basis. We 
will list the coefficients where 

i=l 9 A=l 

The result is 

M) _ 1 p(A) , 1 F (C) _ p (B) _ 3ki 3k 3 7^3 
dl "8 17 + 2^ 17 ^ 28 + ^ 30 167r + 2vr + tt ' 

rf 2 -^10 + ^17 +2^17 +^17 +^28 ~^30 + Zt § 



3«i _ 7^3 _ 47K3 

167T 2-7T 7T 

4 1} = 7^1 + M - k 3 A , 
JX) - 1 F , F (! + 7)«3 

«4 — + ^6 



(3.67) 



4 1} = ("-!) 



4 2vr 
2F 6 + ^' 

7T . 



where 



«3 = k 2 + -«i • (3.68) 

For the calculation of the second term in (|3.59|) we just need ip given by (|3.56|) . 
One then gets: 

(2) n(l - 2-kAq) n 
r n ^T~ 2 + + (3.69) 

Finally the third term can be obtained from the 5x5 two-loop anomalous di- 
mension matrix of the dimensionally regularized scalar operators. The computation 
is described in the last subsection of Appendix A, and leads to the simple result 



.(2)_ o 9 



111 2(n - 2) 

For n = 3 it is -13/2. 



(3.70) 
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Putting everything together, we find for the sub-leading coefficient 



where 



A 2 > = ai + na 2 + — ^ , (3.71) 
n — 2 



19 d 27/ei 3ciKi vrr 3 11 (B) 
8 2 64k 3 8k 3 k 3 I 2 4 4 15 



8^17 -g^!7 "2 F 17 +i F 28 "4^30 +2 F 6) 



(3.72) 



and 



a 2 = — (3.73) 

2 K 3 

a 3 = ---4nA . (3.74) 

Finally we remark that in the n = 3 case it would be nice to know also the 
three-loop coefficient however it would require a lot more effort to compute. 

(2) 

This is built from (among other things) the two-loop coefficient c\ appearing in 

(3) 

the effective action and the three- loop anomalous dimension matrix elements iA. . 
4 Results from Monte-Carlo simulations 

We study here the lattice artifacts of the step scaling function introduced in [21] . It 
is defined as follows. One considers the model confined to a finite (1-dimensional) 
box of extension L with periodic boundary conditions. The dimensionless LWW 
coupling is defined as 

u = Lm(L) , (4.1) 

where m(L) is the mass gap of the theory in finite volume. Next one measures u\ 
defined similarly with doubled box size. In the continuum u\ is a function of uq, 
called the step scaling function u\ = a(2, uq). 

In the lattice formulation one considers an infinitely long strip with L/a sites in 
the spatial direction and tunes the parameter = 1/ Aq so that the measured value 
of uq equals to a prescribed value. (In our case, following [21], we used no = 1.0595.) 
Then with the same (3 one measures 

ui = 2Lm(2L) = S(2, u ,a/L) . (4.2) 

Finally, repeating these steps with finer resolution, one extrapolates to the contin- 
uum limit 

a(2,u )= lim £(2,u ,a/L). (4.3) 

o/i-t-0 
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1.31 



W 1.27 



o 



1.29 




1.25 







0.05 



0.1 



a/L 



Figure 1: Monte Carlo measurements of the step scaling function at uq = 1.0595 for 
two lattice actions in the 0(3) model. The data for the larger artifacts correspond 
to a modified action. The fit contains a and a In a terms. 

The deviation from the continuum limit is denoted below by 



The advantage of this measurement for the purpose of studying lattice artifacts 
is that there is no need to know the box size L or the mass gap m(L) in physical 
units, i.e. one does not need to measure the correlation length in an infinite volume. 
Moreover the continuum 2-dimensional O(n) model is exactly integrable and the 
finite volume mass gap (and hence also the step scaling function) is exactly calculable 
using Bethe Ansatz techniques [22j . 

In the MC simulations we used a modified version of the improved estimator of 
Hasenbusch |24j . This is described in Appendix E. The results of the MC measure- 
ments for S(2, -uo, a/L) for the 0(3) case are shown in Fig. Q]fbr two different lattice 
actions, for the standard one (ST) and an action D(l/3) defined below. One can 
see that the lattice artifacts (cutoff effects) are very nearly linear as function of the 
lattice spacing a both for the case of the ST and D(l/3) action. Although the effects 
are in this case relatively very small, they seem not of the theoretically expected 
form. Note however the encouraging feature that computations with different lattice 
actions are consistent with the same continuum limit, supporting the crucial concept 
of universality, even if both extrapolations would miss the exact continuum limit by 
about 0.002. 



511(2, uq, a/L) = S(2, uq, a/L) — cr(2, uq) . 



(4.4) 
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4.1 The lattice actions 



In the simulations we used a one-parameter lattice action D(a) which includes di- 
agonal interaction: 

A(S) = -(1 - 2a) s ( x ) ■ s i x + «A) - « S ^ ' S ( x + ad ) ' ( 4 - 5 ) 

X,/L« x ,d 

where in the second term the summation is over the two diagonal directions d = (1,1) 
and (—1, 1). 

In the analytic calculations we considered a more general class of lattice actions 
including those corresponding to the kernel with two parameters 

K(p) =p 2 + g ia y + g 2 a 2 [ff . (4.6) 

This includes besides the diagonal interaction also the on-axis second-neighbor in- 
teraction as well. The action D(a) in (|4.5p is a special case with 

9i = ~92 = . (4.7) 

The tree-level Symanzik improved action (SYM) corresponds to g\ = 1/12, g 2 = 0. 
For the in (|2,39p we have 

«i = 91 ~ ' K 2= 92- (4.8) 
The coefficients of the tree-level effective action are given by 

e = -12/ci -16«2 = 1 + 2q, e 4 = -12«i = 1 - 6a . (4.9) 

Results of the numerical simulations for the 0(3) and 0(4) cases using different 
actions are collected in Tables [T]1SJ 

The analytic expression for the lattice artifacts are expressed in terms of the 
inverse lattice coupling, f3 = 1/Ag. To relate the results for different actions it is 
convenient to introduce /3 e g from the two- loop formula (|2,67|) . For the ST action 
one can get this using the relations (|2.68[) . (|2.18p and the results 

ML = 0.2671536 for 0(3), at u = 1.0595 , 

ML = 0.2390557 for 0(4), at u = 1 , (4.10) 

ML = 0.3408255 for 0(8), at u = 1.0595 , 

which can be calculated using TBA techniques [22]. In the range of our simulations 
Peg differs only slightly, by ~ 0.04 from the actual (3. 
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p 


L/a 


Uo 


Ui 


T(2,u ,a/L) 


ST 


1.56923 


5 


1.05947(4) 


1.29375(5) 


1.293798(81) 


0.032588(81) 


1.60476 


6 


1.05952(5) 


1.29018(4) 


1.290148(89) 


0.028938(89) 


1.69850 


10 


1.059484(21) 


1.279805(33) 


1.279830(47) 


0.018620(47) 


1.73020 


12 


1.059495(12) 


1.276719(29) 


1.276727(35) 


0.015517(35) 


1.77900 


16 


1.059492(10) 


1.272430(15) 


1.272442(22) 


0.011233(22) 


1.84603 


24 


1.059492(8) 


1.268074(27) 


1.268087(30) 


0.006877(30) 


1.89295 


32 


1.059498(14) 


1.265932(29) 


1.265935(37) 


0.004725(37) 


2.00553 


64 


1.059516(6) 


1.262938(14) 


1.262912(17) 


0.001703(17) 



Table l: Results of MC simulations for the 0(3) standard action. The last two 
columns are the values extrapolated to uq = 1.0595 and the deviation from the 
continuum value a(2,iio) = 1.261210 [22J. 



According to (|3,55p for other actions one has 

/? eff = /? e s ^ + ^(V-VsT). (4.11) 

Values of ip/{2ir) for various actions are given in Table [6l 

Figs. [2] and 0show the deviations from the continuum limit, ST(2, Uq, a/L), as 
a function of a/L for simulations using different actions in the 0(3) and 0(4) case, 
respectively. 

In the ratio of artifacts for different actions the unknown non-perturbative co- 
efficient Di(A) in (|3.57p . (13.58p cancels. Expressing the ratio in terms of (3^ one 
has 

f ( [ X ° ,a \ =eoU + (p- PST)^ + .-.), (4.12) 

where 

^(" ,21 + ^2*)- ( " 3) 
The values for these constants for various actions are summarized in Tabled 
Figures [H-lSl show the values of ST, (2, uq, a/L)(L/a) 2 vs. the corresponding (3 
values. The 2-parameter fits are the predictions from (|3.57p . (|3.58p . where the overall 
constant and the sub-subleading correction coefficient are fitted, while the value of 
r^ 2 ) is fixed to the known value. 
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Figure 2: <5£(2, u , a/L) is plotted for the 0(3) ST, D(l/3) and D(-l/4) lattice 
actions. 



0.015 



i o.oi 
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•— • 0(4), ST 
• --■ 0(4), D(l/3) 



0.05 



0.1 



a/L 



Figure 3: <5£(2, u , a/L) for the 0(4) ST and D(l/3) actions. 
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Figure 4: <5X(2, uq, a/ L)(L/ a) 2 for 0(3), ST action. The solid line is the fit 
3.1320(/3 3 + c(3 2 - 0.4883^), with c = A 2 ^/(2tt) = -1.1386 fixed. 




Figure 5: 5T,(2,uo,a/L)(L/a) 2 for 0(3), D(l/3) action. The solid line is the fit 
5.4803(/3 3 + c/3 2 - 0.2309/3), with c = r( 2 )/(2vr) = -1.5641 fixed. 
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1.5 



1.6 



_l_ 
1.7 



L.8 



1.9 



Figure 6: <5E(2, u , a/L)(L/a) 2 for 0(3), D(-l/4) action. The solid line is the fit 
1.6566(/3 3 + c/3 2 - 1.7554/3), with c = r( 2 )/(2yr) = -0.1333 fixed. 




Figure 7: <5E(2, u , a/L)(L/a) 2 for 0(4), ST action. The solid line is the fit 0.185 {(3 2 + 
c/3 - 2.170), with c = r( 2 )/(27r) = -0.7295 fixed. 
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p 


L/a 


u 


Ui 


v-l /r-i — / T \ 

E(2,u ,a/L) 


(5S 


2.00460 


10 




1.059521(26) 


1.297390(45) 


1.297357(61) 


0.036147(61) 


2.03654 


12 




1.059578(22) 


1.290733(48) 


1.290609(59) 


0.029399(59) 


2.08453 


16 




1.059485(25) 


1.282422(47) 


1.282446(61) 


0.021236(61) 


2.14870 


24 




1.059544(24) 


1.274298(46) 


1.274228(60) 


0.013018(60) 


2.19293 


32 




1.059526(23) 


1.270176(47) 


1.270135(59) 


0.008925(59) 


2.29813 


64 




1.059515(30) 


1.264571(59) 


1.264547(76) 


0.003337(76) 


Table 2: 


Results of MC simulations for the 0(3) D(l/3) action, 


P 


L/a 




Ui 


E(2,uo,a/L) 




1.60705 




12 


1.059571(17) 


1.270557(37) 


1.270444(46) 


0.009234(46) 


1.65572 




16 


1.059461(26) 


1.268140(33) 


1.268202(53) 


0.006992(53) 


1.72330 




24 


1.059469(24) 


1.265479(35) 


1.265528(52) 


0.004318(52) 


1.83791 




48 


1.059450(20) 


1.262886(101) 


1.262965(106) 


0.001755(106) 



Table 3: Results of MC simulations for the 0(3) D(-l/4) action. 



A further check is provided by the ratio of artifacts. For the 0(3) case the 
ratios for D(l/3)/ST and D(— l/4)/ST actions are shown in Fig. 0where the data 
show a remarkable agreement with the parameter- free prediction in (|4. 12j) . (Note 
however that it is possible that this few percent agreement is due to an accidentally 
very small coefficient of the next, 1//3 2 , term.) 

The corresponding ratio for the 0(4) case is shown in Fig. [TOl Although the 
agreement is poorer in this case, at larger (3 (for L/a = 32, 64) the data seem to 
approach the prediction. 

In [25], [TO] the cut-off effects for the step scaling function were measured for 
0(4) and 0(8) at no = 1.0595 with the ST action. Since our errors for the 0(4) case 
are considerably smaller, we analyse here only the 0(8) data from these papers. 
FigfTT] shows that the MC data of [25j,[10j are also consistent with the analytic 
predictions. (For 0(8) one has er(2, uq) = 1.345757 at uq = 1.0595.) 
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Figure 8: (511(2, u , a/L)(L/a) 2 for 0(4), D(l/3) action. The solid line is the fit 
0.199(/3 2 + cfi + 3.566), with c = r( 2 )/(2vr) = -0.9623 fixed. 
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0(3), D(l/3)/ST 
0(3), D(-l/4)/ST 



I , 1 , 1 , 1 , 1 

1.6 1.7 1.8 1.9 2 

Peff 

Figure 9: Ratio of deviations <5£(2, u , a/L) for the D(l/3)/ST and D(-l/4)/ST 
actions vs. /3|g for the 0(3) model. The solid lines are the predictions 5/3(1 + 
0.224/7^), and 1/2(1 + 0.691//?^). 
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Figure 10: Ratio of deviations <5£(2, uq, a/L) for the D(l/3) and ST actions vs. /3 e g 
for 0(4). The solid line is the prediction 5/3(1 + 0.375//?^). 




30 



(3 


L/a 


u 


Ui 


S(2,M ,a/L) 


8T, 


2.4772 


8 


1.00003(5) 


1.21496(5) 


1.21491(9) 


0.00625(9) 


2.5532 


10 


1.00021(6) 


1.21357(5) 


1.21325(10) 


0.00459(10) 


2.6152 


12 


0.99994(3) 


1.21209(4) 


1.21218(6) 


0.00352(6) 


2.7117 


16 


0.99997(4) 


1.21100(5) 


1.21105(9) 


0.00239(9) 


2.8470 


24 


1.00006(4) 


1.20995(5) 


1.20986(8) 


u.uuizy^o j 


2.9428 


32 


0.99992(3) 


1.20925(7) 


1.20937(9) 


n nnn7i (q) 


3.1726 


64 


1.00002(4) 


1.20887(6) 


1.20884(8) 


0.00018(8) 


4: Results of MC simulations 
= 1 [22]. 


for the 0(4) standard action. 


£7(2, M J — 1.2 


P 


L/a 




Ml 


T,(2,u 0l a/L) 




2.8512 


8 


0.99992(5) 


1.22615(7) 


1.22628(11) 


0.01762(11) 


3.0845 


16 


0.99996(4) 


1.21452(6) 


1.21458(9) 


0.00592(9) 


3.3095 


32 


1.00001(5) 


1.21034(6) 


1.21033(9) 


0.00167(9) 


3.5332 


64 


0.99988(5) 


1.20898(7) 


1.20916(11) 


0.00050(11) 



Table 5: Results of MC simulations for the 0(4) D(l/3) action for mq = 1. 



5 Conclusions 



The main motivation of this work was to explain how the apparently linear artifacts 
found in ref. [2] can be reconciled with standard expectations. We have shown that 
the data, although astonishingly linear as function of the lattice spacing, can equally 
well be described by a more complicated formula we calculated using Symanzik's 
theory of lattice artifacts. 

Although both type of fits describe the data well, we think that by now there is 
no doubt that the conservative description of lattice artifacts based on Symanzik's 
effective action is correct and there is no unexpected new physics behind the phe- 
nomenon (as originally suspected). Our arguments can be summarized as follows. 

1) Since ref. [2] appeared, the exact continuum limit has been calculated [22] 
by bootstrap techniques and using this extra piece of information we see that the 
"linear" fits of Fig. [1] must be bent a little, as seen in Fig. [2j In other words, 
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n 


Action 


^/(2vr) 


A 2 >/(2tt) 


P 


P ~ PST 


3 


ST 


0.37031535 


-1.13861509 


-0.02766903 






D(l/3) 


0.58672731 


-1.56413634 


0.19604560 


0.22371463 




D(-l/4) 


0.26551678 


-0.13328534 


0.66326502 


0.69093405 




SYM 


0.24344345 








4 


ST 


0.49063070 


-0.72946015 


0.25180125 






D(l/3) 


0.79446745 


-0.96233225 


0.62660266 


0.37480141 




D(-l/4) 


0.35356625 


0.22986274 


0.93699525 


0.68519400 




SYM 


0.31658887 








8 


ST 


0.97189212 


-0.58931931 


0.70653685 





Table 6: Values of constants appearing in the description of the cutoff effects for 
different actions. 

continuum extrapolations based on "linear" fits would miss the exact result by about 
0.002. 

2) Lattice artifacts normally decrease very rapidly as a 2 , but here this is par- 
tially compensated by the cubic logarithmic increase of the /3 3 factor for 0(3). This 
increase is further enhanced by the relatively large negative coefficient of the sub- 
leading term. As can be seen in Fig. |4j the logarithmic correction factor increases 
by about an order of magnitude between (3 = 1.6 and (3 = 2.0. This can mimic a 
1/a-like increase in our limited range of (5. 

3) The correctness of the effective action description is further corroborated by 
its parameter- free prediction for the ratios of artifacts, which agrees very well with 
the measured values (see Fig. [9] and Fig. [10] for the less spectacular 0(4) case). 

4) In the large n limit our formula reproduces the results of refs. [10|11|12] . 
After completing this long computation we arrived at the sobering conclusion: 

there are no linear artifacts and Symanzik's theory describes the data well. But 
it was nevertheless useful to go through the steps of the calculation because it 
provided us with the opportunity to learn about Symanzik's theory of artifacts and 
improvement. Similar computations should, in our opinion, accompany precision 
lattice studies of QCD in order to control and better estimate systematic errors 
arising from lattice artifacts for extrapolations to the continuum limit. 
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Appendix A: Anomalous dimensions of the dimension 4 operators 
A.l Scalar operators 

The one-loop mixing matrix for dimension four invariant Lorentz scalars has been 
first computed by Brezin et al. [18] using the basis (l23H|) 1. It is given by (|2.26l) 



withB 



2nwij 



4 - 2n 


-4 





n 
2 


-2 


2 - 2n 





n+2 
4 


4 


-16 


2 - n 


-2 


-4 


16 





4 - n 


1 


-4 


l 

4 









4(n- 1 



(A.l) 



/ 



In Fig. [T2] we show the diagrams which are needed for the renormalization of oper- 
ators 0\ and O2 to 1-loop order. 




>-8- <> -8_«= 



Figure 12: 1-loop contributions to the renormalization of operators 0\ and Oi- 
Full circles represent the (4 or 6-leg part) of the operators, full squares stand for the 
4-point interaction vertex of the model. 



We will now simplify the mixing problem using some operator identities. The 
identities can be derived by making in (|2.3j) the infinitesimal change of variables 



6^ = SujQ 1 , (A.2) 



3 where, for simplicity, all operators can be taken at zero momentum (space integrals) 
4 (and we have checked their result) 
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where 5uj is an infinitesimal parameter and Q l is an arbitrary local expression of the 
fields and sources. The corresponding change in the action is 

5A = -Su J d D x (on* - an 1 + \ Q i (A.3) 

and since the measure (T>7r) is invariant in dimensional regularization under any 
local change of variables we have 

Z[I] = J (Vtt) e- A - 5A = Z[I] - j (Prr) SAe~ A , (A.4) 

leading to the operator identity 5 A = 0. Thus we have 

□tt* - + —A Q { ^0, (A.5) 

where ~ means operator identity for the corresponding space integrals. 

We now consider the infinitesimal transformations corresponding to respectively 

Q i = {go (d»S ■ d^S) ,g 2 (I-S + a) ^\g I l + g 2 av\ g 2 U^} (A.6) 

and find the identities 

(i ■ s) (d^s-d^^ad^s-d^s, 

(/■S)W, 

i 2 + i-os+(i-s) (d^s ■ d^s) « a 2 , 

I -OS- (d^S ■ d^S) 2 + DS ■ OS « -a^S ■ d^S) , 

respectively. Here we have introduced the operator 

a = a + d tx S-d lx S. (A.8) 

We see that insertions of the apparently O(n) non- invariant operators O4 and O5 
are equivalent to inserting the manifestly invariant ones in (|A.7p . There are no 
further identities independent of the ones in (IA.7jl . This can be shown by using the 
continuum analog of the lattice considerations of ref. [5] . 

It will be convenient to use a new basis of operators. We keep the "hard" 
operators 0\ and O2 but instead of the rest we first introduce W, U and V, where 

AW = ad tx S-d li S, 

8U = a 2 , (A.9) 
2V = OS ■ OS — (d M S • d^S) 2 + a 2 . 
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This basis change corresponds to 

e> 3 = V - 4Z7 + 40i , 

O i = 2W-40 1 , (A.10) 
O b = U -W + O x . 

Further we introduce the combinations 

U 3 = W-2U + \v, U± = U--^V, U 5 = V. (A.ll) 

2 4n 

Now the operator identities can be rewritten as 

c7 4 » Inv 6 , (a. 12) 

8 

& 2 

r afe = S a 5 b - -5 ab . (A. 13) 

n 

Inserting the above identities into the generating functional (12.191) we can derive 
useful identities for the correlation functions of the operators Uj ,j = 3,4, 5. These 
are best formulated in Fourier space. We define, as usual, 



where 



J d D xe iipiXl+ - +PrXr) £^- ar (xi, ... ,x r ) 



(2Tr) D 5(p 1 + ---+ Pr )g^ n - a r(p 1 ,...,p. 



(A.14) 



for correlation functions with and without operator insertions. With this notation 
we have 

%^( P l, . . . ,p r ) = f OTpI) d"'( Pl , ..., Pr ). (A.15) 



\k=l 



From this formula it is clear that the operator U$ is renormalized multiplicatively 
with renormalization constant 

Z V =^ = l- (2 - n)g2 +... (A.16) 
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For t/4 we have 



u 4 









G* iar > a3 - ar (pi+P2;P3,- 


■ ,Pr) + 


[(;)-] 









2 2 

+ ^ g° in < a2 - ar ( Pl ;p2, ■ ■ ■ ,Pr) + [r-l] perms 



perms 



in 



nn\a\...a r 



(0;pi,...,p r ), 



(A.17) 



where 

gf ;ai - a -(p;w,...,p r ) (A.18) 

is the Fourier space correlation function of the local operator r ab at momentum p. 
Actually, (|A.17p is valid for r > 2 only. For r = 2 we have 



^; 2 b, -p) = ^<r^> + zMgpwfa _ p) 



m 2 ^ 



m 



(A.19) 



+ -^<? r a2n;ai (-p;p) + — er ;aia2 (0;p, -p) • 



Since m (|AT7l) (and in ([ATL9]) ) all terms require the same overall renormalization 
constant, the operator [^4 renormalizes multiplicatively with renormalization con- 
stant 

ZZ T _ (3 - 2n)g 2 



u * Z 2 
where we have used the result 



1 



2vre 



+ ..., 



(A.20) 



Z T = l + ^- + ... (A.21) 
lire 

Finally for U5 we find 

G^- ar (pi,---,Pr) = (2TT) D 5( Pl + P2 )g*g a ^ Q -1 perms 

+ {2n) D 5{ Pl )5 a ^ n m 2 g 2 g a ^(p 2 , . . . ,p r ) + [r - 1} perms 



(A.22) 



for r > 2 and 



G^ 2 (P, ~P) = 9& aia2 + 2m 2 g 2 5 a ^ n (a)5(p)(27T) D 



(A.23) 
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for r = 2. Again, U$ renormalizes multiplicatively with 

Z , (3 - n)g 2 



% 5 - Z2 - ' 



2vre 



+ ... 



(A.24) 



We can check the results (lAlBl) . (|A~20l) and (|AT24|) by using the mixing matrix 
(lA.lh , We define the linear operator 

C(Oi) = 2-KWijOj (A.25) 

and then find 

C(U 3 ) = (2 - n)U 3 , C{Ui) = (3 - 2n)U 4 and C(C7 5 ) = (3 - n)F 5 . (A.26) 

There are two other linear combinations Ui,i = 1,2 diagonalizing the one- loop 
mixing matrix (jA.lj) with 



C(C/i) = 2C/i and C(J7 2 ) = (4 - 2n)C/ 2 

The relations between the bases are given by 

Ui — SijOj , 



(A.27) 



(A.28) 



with the matrix S given by 
/ n+4 



S 



2n-l 

5 -2 



1 



(n-2)(n+l) (n-2)(ra+l) 2(n-2)(n+l) \ 
2n(2n-l) 2n(2n-l) 2n-l 



1 

"n-2 

1 
2 

4n 



2n-5 
n-2 

1 

2 

n-1 
2n 



4 



n-1 











(A.29) 



For later use we also write down the inverse of this matrix 



s- 1 



( _JL_ 


1 


n 2 -2n-4 


2(n+4) 


(n 2 -n-4) 


n-1 


n-1 


n(n-2) 


2n-l 


2n(n-l) 


1 


n 


n 2 — 4n— 4 


5n+2 


(n 2 -3n-2) 


2(r»-l) 


2(71-1) 


2n(n-2) 


2n-l 


4n(n— 1) 


4 


4_ 


4(n 2 -2n-4) 


36 


(n 2 -9) 


n-1 


n— 1 


n(n— 2) 


2n-l 


n(n— 1) 


4 


4 


2(n 2 -2n-8) 


36 


n 2 -9 


n-1 


n— 1 


n(n-2) 


2n-l 


n(n— 1) 


V n^T 


1 


4 


9 


(n-9) 


n-1 


n{n-2) 


2n-l 


4n(n— 1) 



(A.30) 



37 



A. 2 Tensor operators 

Here we discuss some properties of the four-index symmetric tensor operators denned 
in (gjEH) and (|23Tj) . 

Let X pupa be a tensor operator, totally symmetric in its indices fi, u, p, a. We 
can construct the corresponding totally symmetric, traceless tensor operator X pupa 
by subtracting trace terms as follows. 

X pvpiJ = X pvpa 4-4 {sV u/ £/ 3fT 9po£,p.v ~i~ 9pp£,va ~\~ Qva^pp 

r (A.31) 

"I" 9p.u£,vp + 9vp£,p<j j" "I - ~ 2)(Z) + 4) \p^ v ^P a 9pp9ua + 9pa9up^ > 

where 

Cpo- = ^aopcr, and £ = X ao ,pp . (A. 32) 

In the effective Lagrangian there are operators of the form X^=i X pppp > which can 
be rewritten, using (|A.31|) as 

D D 3 

^2 X ww = ^2 ^ww + 4^17 X ""/3/3 • (A.33) 
If we apply this procedure to the operator t pvpa in (|2.30p we get 

J2S-df l S = A + -—0 3 , (A.34) 

1 

where ^4 is defined in (|2.32p . The traceless symmetric tensor operator i^p^ can mix 
under renormalization with other traceless symmetric tensor operators only (which 
must have the same values for other quantum numbers). All tensor components 
renormalize the same way and using this property we can simplify the renormaliza- 
tion problem by considering the tensor component, where 

V+ = Vi + iV 2 (A.35) 

for any vector. Great simplification occurs in PT calculations for this tensor compo- 
nent since g ++ = and many diagrams vanish, and it is thus easier to renormalize 

A = t ++++ = d\S ■ d\S (A.36) 

instead of A itself. (Note that we are considering operators at zero momentum, i.e. 
their space integrals.) 
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tfiupa can mix under renormalization with k^ vpa given in (|2.3ip . For this oper- 
ator we find 

D 8 

(d^S ■ d^S) 2 = B + — {ft + 20 2 ) , (A.37) 

where B is defined in (|2.33p . For convenience we define 

B = k ++++ = (d+S ■ d+S) 2 . (A.38) 

There are no other spin-four, dimension-four O(n) invariant operators, there- 
fore we have a 2 x 2 renormalization problem here. By calculating the matrix 
elements of A, B we find that the mixing matrix is diagonal, 



g 2 (n - 2) /0 
2ire 



% = <% + \„ ' I I + • • • , (A.39) 



,0 1, 
for the operators 

U 6 = \a, U 7 = \b . (A.40) 

% % 

A. 3 2-loop Mixing matrix elements 

(2) 

For the computation of the mixing coefficients zA , we worked in the Brezin et al 
basis (|2.28|) and computed the corresponding pij coefficients of the g 4 /e terms in 
(|2.26p . From the relation between the bases (|A,29|) the v\- are then given by 



v 



(2) (n-2) 
« (2tt) 2 



ko + S ikPkl (S-% . (A.41) 



For the computation of 7j n we only require v\- for i,j = 1, 2. For this purpose we 
only need the for i = 1,2 because using the fact that Kij = for i = 3,4, 5 and 
j = 1, 2 we have 

5 2 

Y^SikipS' 1 )^ =J2 7 ik(pS~ 1 )kj , {1,2}, (A.42) 

k=l k=l 

with 

5 

%k = Sik - Yj Sir<JrsS s k , £, fee {1,2}, (A. 43) 

r,s=3 

where a is the inverse of the lower 3x3 diagonal block of S (i.e. X^s=3 a rsS s t = $rt)- 
For the computation of p^j for j = 3, 4, 5 we need to compute the divergent parts 
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of insertions of Oi in the 2-point pion correlation functions (with just up to one 
interaction vertex), and for the computation of pij for 1,2 we need the analogous 
computation for the insertions of Oi in the 4-point pion correlation functions (with 
up to two interaction vertices). 

The computation is again too lengthy to present the details here; we just list the 
relevant Feynman diagrams in Fig. 1131 Note that the 1-loop contribution correspond- 
ing to the diagram b in Fig. [12] is only a wave function renormalization. Analogous 
contributions (with wave function renormalization on one or two external legs) are 
not shown among the diagrams of Fig. [T3l Diagram d in Fig. [12] contains an internal 
contraction line and thus vanishes in the limit m — > 0. Diagrams containing similar 
internal contractions are not shown in Fig. [13] either. 

We checked that the 1/e 2 divergences were as expected from the RG consider- 
ations, and from the 1/e divergences we finally obtained 



> I 7n n n — 3 
(2ir) 2 Plj = \ _-H ) -9n + 19,-,^-,-n + l 

> i 13n _ 3n — 13 n + 2 n — 1 n(n — 1) 
(2,? P23 = {-— +7 t —^ -—,--L-J. 



(A.44) 



This yields 



(2) = _3(2n»-5n + 5) 
v ; 11 2(n-l) 

fn\ nin - 

(2ir) 2 v, 



2 (2) _ n(n - 4) 
21 ~ (n - 1) 



(A.45) 



Appendix B: Lattice perturbation theory 
B.l Perturbation theory with periodic be 

We consider a volume V = N 2 with periodic boundary conditions in each direction 

S(x + Np) = S(x) , n = 1, 2 . (B.l) 

Let O be an O(n) invariant observable. For Ao — > zero modes appear due to the 
O(n) invariance. To avoid this Hasenfratz [23J used the Fadeev-Popov trick: 

(O) = \ J d»r J [dS]e- A 6 (r-^ S(x)j 0[S] . (B.2) 

Because of the O(n) invariance the inner integral is independent of the direction of 
r. Take r = \r\n$ , no = (0, 0, . . . , 0, 1) and set 

S(x) = (X TT(x),a(x)) (B.3) 
a(x) = y/l - A2vr(x)2 (B.4) 
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where w(x) has n — 1 components 7r J (x) , j = 1, . . . ,n — 1 . Then 



-— / 



5 ^tt(z)^ Oe^* , 



(B.5) 



with 



4ff = 4 + ^lna(z)-(n-l)ln(— ^ff(a;)] . 

X \ X / 



(B.6) 



Then when one expands for Ao — > no zero modes occur because of the delta 
function. 

In the following we restrict results to the infinite volume limit. 
B.2 1-loop connected 4-point function 

The leading orders of the connected 4-point function have been given in subsect. 2.2.2. 
Here we give the result for C2 and consider the coefficients of powers of (n — 1) sep- 
arately: 



c 2 = 5>-iycT. 



(B.7) 



r=0 



~[2l 

First for C\ 1 we obtain simply: 



C2 ] (pi,qy,P2,q2) = 1 



2 K K K K 



K pi+qi F(p 1 + qi ) . 



(B.8) 
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Next for we get: 



IK K K K 



[ K pi+ P 2 F (Pi +P2)+ K pi+q2 F{p 1 + q 2 )] 



- l - [F(pi) + F( qi )) C?\p u qv,P2,<!2) ~ \ [F(P2) + F(q 2 )} cf\p } .qy. p 2 . q 2 I 
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2 K P2 Kq 2 
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a; 
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+ 77 



2K P1 K 91 
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K K K K 



Pi iv 9l J 



^(P2) + F(q 2 ) 



K 



pi 



K, 



92 J 



K K K K 



-if, 



Pi+qi 



k 



Pi 



iKpi 
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k, 



91 



K„ 



+ 



k 



■p 2 



+ 



"92 



+ 



pi+qi 



K K K K 

i\p 1 i\q 1 i\p 2 i\q 2 



qi ^P2 
\Ji(pi,Qi) + Ji(P2,q 2 )] 



K, 



92 



K K K K 
1 ^pi iy qi 1 ^P2 1 ^q2 



\Z A {pi,p 2 ,qi) + Z 4 (p 1 ,q 2 ,q 1 )] 



-K, 



pi+qi 



Ti(pi,qi) T 1 (p 2 ,q 2 ) 



K K 

i\p 2 i\q 2 



Kpi Kqi J 



(B.9) 



where 



Z 2 (p,q) = [ 



p-q 



K 



p+q 



_ Kk Kk-p K k K p KpKk-p 
Jl(p,q) = Z 2 (-p-q,p) = Z 2 (-p-q,q) 
K k K„ 



Ik 



Kk-pKk+q K p+q Kk-p Kp +q Kk+q_ 



(B.10) 



(B.ll) 



-7 i \ f \ K k-pK k 
Jk Y K k^k-p 

Ti(p,q)= f 
Jk 



k+r 



K p K r 



Kp-\-q-\- r Kq 



-Kk-Kp+q KpJ r qKk—p—q_ 



(B.12) 



k [KkKk-pKk+q K k KpK q KpK k _pKp +q KqKp+qKk+q 

(B.13) 
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Finally for C 2 we get: 
1 



C l 2 (Pl,9i;P2,<?2) 



K 

P2- tl 92 L Iy Pl 



F(pi) + F(gi) 



A', 



+ 



xv pi xv <ji 



^(P2) + 



P2 



K, 



92 J 
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+ 
+ 



1 1 

+ 



K K K K 

1 1 



^Pl+P2 , ^Pl+92 
If 2 + 7^2 

P1+P2 Pi +92 



1 



K K K 



kq2 _|_ ^P2 



l '/2 



K, 



P2 



+ 



1 



k qi k Pl 



l 91 



K, 



pi 



K K K K 



[^pi "I - "I - kp2 ~t~ ^92] 



K K K K 



[Kpi+P2 + Api+92] 



V P1 



^<2l _|_ ^P2 _|_ ^92 



A' K K K 

_ V P1 J1 9l J1 P2 - fi 92_ 



1 



^P2 ^92 . Api +P2 



1 



{ ^(Pl, P2) + Z 2 (qi, Q2)} + 



K 



K Pl K qi 



K, 



{Z 2 (p2,Pi) + Z 2 {q 2 ,q{)} + 



{Z 2 {pi, Q2) + Z 2 {qi,P2)} 

pi+92 

1 



K r , 



+ 



K K K K 



P1+P2 -"Pl+92 

^pi+p2 {Ji(pi,P2) + Ji{qi,Q2)} 



{Z 2 {p 2 ,qi) + Z 2 (q 2l pi)} 



+ Api+92 {Ji(pi,Q2) +Ji(qi,P2)} 
1 



+ 



Kpi Kqi Kp2 Kq2 



Zi{pi,qi,P2) + Z A (p 1 ,q 1 ,q 2 ) + Z A (p 1} p 2 ,q 2 ) + ^4(^1,^2,^2) 



1 



K K 
ly P2 IY q2 



\T 2 {pi,qi,p 2 ) +T 2 (p 1 ,q 1 ,q 2 )] 



K K 



+ B(pi,qi,p 2 ) + B(pi,qi,q 2 ) + B(p 1 ,q 2 ,p 2 ) 



[T 2 (p 2 ,q 2 ,pi) +T 2 {p 2 ,q 2 ,qi)] 



(B.14) 



where 
T 2 (p,q,r) 



K, 



k—p—r 



K, 



p+r 



KkK k -pK k+q K k KpK q 



K r 



K, 



p+q+r 



KpKk-pKp +q KqK p+q Kk+ q 



, (B.15) 



B(p,q,r) 



KkKk-pKk+qKk-p-r KkKpK q K p+r 

1 



KpK k _pKpj r qK r K n Kr,+ n Kk+nK,, 



qJ^-p+qJ^k+q^p+q+r 



p+r J^r^p+q+r K- k—p—r 



(B.16) 
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B.3 Lattice integrals 

The expansion in the cutoff of the functions F(p),k p and others appearing in Ci 
(which will be considered in the next subsection) involve various lattice integrals 
that we will first define here. They depend on the specific lattice action and hence 
on the associated function K p (12.371) . We first define 

A o S^ = ~J^[K^(s)-K^(0)} , (B.17) 

1 f 1 

Atf^xp + -A 2 s^ Xp = J ^- [K^ P x(s) - K^pxtO)] , (B.18) 

where ^...^(s) =8^... d pn K Sl 5^ Xp = S^S^xhp and 

Sp,vX P = Sp,v8\p + 5p\&vp + 8p, p 6\ u . (B.19) 

Note 



K^(0) = 25^, 
Kp.vp\(fy = 8a 2 [3ki5^ap + «2V*p] . 



(B.20) 



The values of A{ and other constants are given in Table [7] for various actions. 
Introducing #b(/c), which restricts the momenta to the Brillouin zone: 

2 

--\\9{TT-a\k, L \) (B.21) 
p=i 



we define 



fMfc) i } 

\ K k A; 2 (a 2 fc 2 + 1) J 



Cl = -4« r^- —±— , (B.22) 



2 /■ fgfl(fc) 1 2a 2 r(AQ 
6 " ./ool *2 (fc 2 ) 2 + [fc 2 ] 2 (a 2 fc 2 + 1) ' ' (L> - >> 



where we have introduced the shorthand 

d 2 k 



(27T) 



2 ' 



(B.24) 



Also introduce 

c 2 = 2vrF 6 - — Ki . (B.25) 



Further constants appearing are 



F M = a~ 2 



g B (fc){jjr w (fc)-.g w (0)} 2 (12«i + 16k 2 ) 
K 2 a £; 2 (a 2 A; 2 + l) 



(B.26) 
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Fi^fiuXp — a 2 



I DO 



[k 2 ] 2 (a 2 k 2 + l) 



and 



-F [A) 5 , +F {B) s s 



(B.27) 



Jk K k 



V 28pvXp — a 



— 17 °p,v\p ~r ^17 u fiu (J Xp ~T r \i ftpj/Ap • 

d B {k)Ky{k)K v {k)K x {k)K p (k) lQk p k v k x k 



(B.28) 



^ 4 



^29 = a~ 2 



6 B (k)K p (k)K p (k) 



I CO 



4 4a 2 r(/c) 

[P]2 " [Pp(^2p + i) 



- ^a/3 7 r(0) ^2i5(° 2fc2 7 + fj \^pk v k x k p + 3 perms] - 2k T kfj,k u kxk p } 

— ^28 u pvXp i r 28 b pvXp > 

(B.29) 



(B.30) 



F^OpvXp — a 



-2 



9 B (k)K„(k)K u (k) 



{K Xp (k) - K Xp (0)} 



2^a/3Ap(0) 



k a kfjk p k v 



[k 2 f{a 2 k 2 + 1) 



(B.31) 



30 °p^Ap 



+ F^buu&Xo + F, 



30 °py°Xp 



30 *p^Ap • 



Note that there are various relations among the integrals defined above. We 
have not attempted to determine all of them but note the following identities 



3K 



F 



2s 



(A) 
30 



29 



3K 



(C) 



30 



30 



8 ^-^ ( 5 B) + ^i + ^(7-4 Cl ) K2 



16vr 



2tt 



(.4) 



2S 



47T 



1 2 

8Ai + -F^f + - (ci - 1) «i = , 

3 7T 

8A 2 + F^ ) + - (c x - 1) k 2 = . 

7T 



(B.32) 

(B.33) 
(B.34) 
(B.35) 



These appeared as consistency relations in the course of our calculation and also 
serve as useful checks on the numerical evaluation, and one can check are satisfied 
by the values given in Table [71 
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B.4 Integral expansions 

Here we give, without derivation, the expansion of the lattice integrals to O(o 2 ). 



HP) = ^ \L{ap) + ci] 

i 2 2 1 



5 p 4 
+ k 2 j L(ap) + -Kipp + c 2 



A; p = ^loP 2 + « 2 [Aip 4 + ^ 2 (/) 2 ] + 0(a 4 ) 
The constants in (|3.3[) are thus given by 



+ 0(a 4 ), (B.36) 



(B.37) 



c 3 = 47r^i + m I - - c\ - 8ttA 
c 4 = c 2 + 4^2 - k 2 {ci + 87rAo) • 



(B.38) 
(B.39) 



~k7 + 4 Fl0 



967rp 2 



K^ypxityq^ -3(p + g)„ (2p p p A + <5 pA p 2 {^(ap) - 1}) 



+ 2 {L(ap) + ci} {(2p v + 3q u )p p p\ + q v q P {2p\ + q\)} 
To consider Z4 we first define the function Z5 through: 
Za(p, Q, r) = Z 5 (p + q, -p, r) . 



(B.40) 



(B.41) 
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Its expansion is given by 

ZbijP, Q, r) A (p + q + r) 2 + q p r u i 4F 13f , u (p) - 8^ 



(c i ~ 1) 
2vr 



F(p) 



+ [{r 2 +pr)(q 2 + pq) + {pr)(pq)] 
- a 2 [A, { (p + qf + {p + rf-p A } + A 2 { [(p + g) 2 ] 2 + [{p + r) 2 ] 2 
+ 24^Ap(0)|4 [prg M + pgr M ] p^ A Pp 
+ [6 {r p r u pq + q^q v pr) + 2 (r 2 ^ + <? 2 r^) p„] p A p p 
+ 2 [r^r ! ,r A (g 2 + 2pq) + q p qvqx{r 2 + 2pr)] p p 

+ r p r u rxr p (q 2 + pq) + q p q u q\q P (r 2 + P r )}^^2 i L ( a P) + c i] 

+ ^ | ^19 M i/(p) + ^22^(p) j + ^ [r>7V9A + q^rx] F 18pu x(p) 

+ ^ [^r^rA^p + q^quQxrp] [2K pu x T (0)F 13rp (p) + a 2 Fi 5 ^ A p] , 

+ -^r^qxqpFi^xpip) . 
The various expressions appearing here are given by 



1 



PpPx + 2 5 P xp2 ^ ~ L ( a P)} 



47rp 2 

Fi6p V \ p (p) = a 2 [25p U 5x p F 10 + F 17pu x p ] 

+ i- [5^V ff (0) + <5 Ap *W(0)] 



[L(op) + ci - 1] - -<5 r(T [L(op) - 1] 



r ■ ' - • - 2 

^18^aO) = ^uPpKx P ra(0)F 13T(T (p) 

+ ^pr(O) [pA^13pr(p) + 2V F 13Ar(p) + PrF 13p x(p)} 
+ a 2 p r ^^Ar^lO + F 15 p uX r + i^i/Ar] , 



(B.42) 
(B.43) 



(B.44) 



(B.45) 



Fl9ixv{p) = [K^pXr(0)P,y + K vp x T {0)p P } P\F 13pr (p) 

+ a 2 pxP P [2S p xSu P F w + Fi 7p xu P ] , 



(B.46) 



7 22 P v(p) = g^aApr(O) [5 U aF 2 4 p X P T (p) + (^a^z/Apr fa)] 



-4F 25/ll ,(p) + -pxP p F26pv\ P 



(B.47) 
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with 



F24a\pr(p) = J [PpPrF 13aX (p) ~ PaP\F 13 pr(p)] + F 27aX pr(p) , 



(B 



24 



24 
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+ PmPi/ <j 2L(ap) - — + ^ (Pp + " 4 



(p 



,2\2 



^ [p 2 <V + 6p,*Pi/] , 



(B 



F26pu\ P = F 2 spu\p - 8\ P 8 f j,vF 2 9 - F 30 ^ u \ p - [5^\5 vp + 8 W 8 V >\ [2F 6 + Fio] 
+ Fi 5fiu \ p - -Fi 7fl \ up — -Fi7n pi/ \ ' (B 



F27p,up\(p) — 



96vr 



L i a P) ~ ^ ) SpupXP 2 



^L(ap) - ^ {pp,Pu8 p x + 5 perms} - s PtlP,/ ^ pPx . (B 



For the scaling part of Z4 this means (p + g + r + s = 0) 

1 



^4(p,9,r) = -A)(a + r) ' + 



vr(p + q) 2 



[(pr)(qs) - (ps){qr)} 



+ 



1 



2vr(p + g) 2 
+ 0(a 2 ). 



{L(a(p + q)) + ci} [(pg)(rs) - (pr)(qs) + (ps)(gr)] 



For the first triangle diagram integral we have: 

1 r. , , . , 



Ti(p,q) 



with (here r = — p — q) 



2irp 2 q 2 (p + q) 2 



(B 



[t (p,g)+a 2 ti(p,(?) + O(a 4 )] , (B 



*o(P, q) = -qrL(ap) - prL(aq) - pqL(ar) + - (p +q + r )ci , 



(B 
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hip, q) 



1 

"2 C1 



p 4 (g 2 + r 2 ) g 4 (p 2 + r 2 ) r 4 (p 2 + g 2 ) 

H o 1 n J> «l 



p* g z 
+ 2(p 2 g 2 + g 2 r 2 + r 2 p 2 )K 2 



-K\ + K2 ) [p 2 g 2 £(ar) + r 2 q 2 L(ap) + r 2 p 2 L(ag)] 



+ 



«1 



(p2)2 4 



{p 2 q 2 V5{p,q) + pVV 5 (p,r)} 



(B.55) 



+ 



7 4 3 



(g 2 ) 2 4 
r 4 3 



{q 2 p 2 V 5 (q,p) + g 2 r 2 V 5 (g,r)} 



-2\2 



- } {r 2 p 2 V 5 (r,p) + r 2 g 2 V 5 (r,g)} 



Here 
V 5 (p,q) 



{(p 2 f + (q 2 r}L(a(p + q ))-(q 2 ) 2 L(aq) 



b 2 ? 2 -2(pg)2] 

{ (p 2 ) 2 + 2p 2 g 2 - 4(pg) 2 } L(ap) + 4p 2 g 2 - 2pq(p 2 + q 2 + 4pg) 



(B.56) 



Note for momenta such that p 2 q 2 = 2(pq) 2 one has g 4 /(g 2 ) 2 -3/4 = -p 4 /(p 2 ) 2 + 3/4 
so that t\ isn't singular in this case H. 

For the second triangle diagram integral we first define 



T 2 (p, q, r) = T 6 (p, q, -p - r) . 



(B.57) 



Then 



Te(p, q, r) ~ F ^ + ^ + r M T 6/x (p, g) + ^^rJT^p, q) 



K 



p+q 



1 r i r 

+ 6 r ^^Ap(0)r A |^— {L(a(p + g)) - L{aq) + L(ap) +ci} 

-(p^g) +-r p Ti(p,g)}, 



(B.58) 



with 

Tenv(p,q) ~ 25,j, v Ti(p,q) 



+ jp -> g| - G 4p a(p,c 



5 indeed we have the identity j-%t 



H(P + 9) 4 + (p-g) 4 -2p 4 -2g 4 }] 



(p 2 ) 2 1 (<? 2 ) 2 



(p 2 <J 2 ) 2 



x [§pV + {pqf 
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where 

G^ip, g) 



8vr(p + g) 2 {P + ^ ~ {p2 + q2 + Apq W a ( p + «)) 

- ^V 2 (p,g) - ^V 2 (g,j>) + ^(p 2 -pq)V 3 (p, q) + ^(q 2 - pg)V 3 (g,j>)} 
|_2 P 2 + 2pqL{a(p + g)) + V 2 (p, g) - (p 2 - 2pg)V 3 (p, g)} 
{-2g 2 + 2pqL(a(p + g)) + V 2 (g,p) - (g 2 - 2pg)V 3 (g,p)} 



+ 



+ 



q^qv 



- (Pfiqu + g^Pi/) <! 2 + -V 3 (p, g) + ^V 3 (g,f>) 



where 



V 2 (p, g) = (P 2 + g 2 )^(«(p + g)) - (p 2 + 2pq)L(ap) - q 2 L(aq) , 



V 3 (p, g) 



(B.60) 



(B.61) 



{4(pg) 2 + 2p 2 (pg) -p 2 g 2 - (q z f}L{a(p + g)) 



,2\2-| 



[p 2 g 2 — 4(pg) 2 ] 
+ {p 2 g 2 - 2p 2 (pg) - A(pq) 2 }L(ap) + (g 2 ) 2 L(ag) 

+ 4(pg) 2 + 2g 2 (pg)-2p 2 g 2 ^ 



Finally 



(B.62) 



T^(p,q) 



where 



p+q 



[G 6fM (p, g) - 4G 3At (p, g) + Z 2 » - Z^(q)} 



Z 2 ^{p) = jP^Fw + o 



2 



F(p) 



GV(p, g) = -i- [L(a(p + g)) - L(ag)] - (p ~ g)} 

GV(p, g) = GV a/1 (p, g) + G 7 6 M (p, g) + GV c/1 (p, g) 
-\(jP-q)» {F W + 4F 6 -F 29 } , 



(B.63) 



(B.64) 

(B.65) 
(B.66) 
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with (t = p + g)): 



G 7a p(p,q) 



1 



+ 



24vr 
PvPxPp 



p* 
(p <-> g) 

1 



^s^Paj jv 3 (p,g) - L(ap) + || 



(B.67) 



G 7bfi (p,q) = —K vXpT (0) -PM^xprfaq) 

+ ( s ^A P Pr + 4 perms) ^L(ap) - ^ 

+ ^2 (<WWr + 9 Perms) + 4 ^ 2)2 P^WPt ^L(ap) - ^ 
- (jp <-> g)} 



IT 

6, 

(B.68) 



GVquO, ?) = ^^Apr(O) t 2 {M^ Apr (p, -i) - PpG 5uX pr(P, ~t)} 

+ t -^ L ^{L(at) + L(ap)-L(a q )} 
+ — rj: {s^a P tPm + 4 Perms} j-L(ap) + ^ j 
j {SfMvPxPpPr + 9 perms} + J 1 2 {^aP p Pt + 5 perms} 



247rp : 



32vr 
- (p <-> g) 



167rp 2 

1 P^PuPXPpPr 
12tt (p 2 ) 2 



(B.69) 



where 



M^px T {p, q) = 47r (J~ +q y {~^ i s ^pXPr + 4 perms] V 3 (p, g) 

VpAr(p)^ 

- (p <-> g)| , 



+ Wr(p)Vio(p, Q) + P ^ 2 P y Pr {V 3 (P, g) + L(o(p + g))} 



(B.70) 
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where Ip U pXr(p) is the totally symmetric traceless tensor: 



PpPvPpPXPr 
(p2)2 



8p 



2 [SpvPpPxPr + 9 perms] 



+ |g [s^pXPr + 4 perms] 



and 



V w(P,q) = ^ p2q2 _ 4 ^ 2 j [(p2 g 2)2 _ I2p2q2(pqy + 16(pg) 4 ] I 

(g 2 ) 2 [(pY ) 2 - 2p 2 ((7 2 ) 2 M - 12pV Cp^) 2 + 8g 2 ( M ) 3 + 16( M ) 4 ] L(a<?) 
-p 2 [(p 2 ) 3 (^ 2 ) 2 + 2{p 2 q 2 ) 2 pq - 8{p 2 ) 2 q 2 {pq) 2 - 2Ap 2 q 2 {pqf 

+ 16p 2 (pq) 4 + 32(pq) 5 ^L(ap) 

+ (P + q? [CpV) V " q 2 } + 2p 2 (q 2 fpq + 8p 2 q 2 (pq) 2 {q 2 - p 2 } 
- 8q 2 (pq) 3 (p 2 + q 2 ) + 16p 2 (pg) 4 l L(a(p + q)) 



-2pq 



(P 2 ) 3 (q 2 f +P 2 (q 2 ) 2 pq{2q 2 + 5p 2 } + p 2 q 2 (pq) 2 {6q 2 - 8p 2 } 



2„2, 



- 8q 2 (pqf{q 2 + 5p 2 } - 8(pg) 4 {g 2 - 2p 2 } + 48( M ) 5 ] } , 

and 

G^ U px{p,q) 



1 



-^s fMUpX pqL(a(p + q)) 



8irpq(p + q) 2 
PuPuPpPx 1 \ 1W s / 

-gV^jw?)-! 



qp,qvq P qx 



(P 2 f 



((?2) 2 

+ {-Hju^aCp) + V4(p, (?) 



where H pup \(p) is the totally symmetric traceless tensor: 
P/xPuPpPx 



and 

V 4 (p,9) = - 



(p2) 



2^2 



(p + q) {pq(p + ? ) - p q ) L(a{p + q)) 



(B.71) 



(B.72) 



(B.73) 



(B.74) 



bV " 2(pq) 2 ] V- 

+ P 2 {p 2 q 2 -pq(p 2 + 2pq)} L{ap) + g 2 {p 2 g 2 - pq(q 2 + 2pg)} L(ag) 

+ 2pq{2p 2 q 2 - pq(p 2 + q 2 + 4pg)}l . (B.75) 
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Writing 



we have for the scaling part 



t 0) + a 2 7t 2) + ... 



Tf\ P , q ,r)=T>( P , q ,r) + ^\^ + 



+ 



{p + q + r)' 

4-7T I p 2 q 2 p 2 (p + q) 2 q 2 {p + q) 2 
where T' 2 is the corresponding continuum triangle integral defined by 

(k — p — r) 2 

(p + q + r) 2 



T^p,q,r) 



k 2 {k-p) 2 {k + q) 2 
(p + r) 2 r 2 
k 2 p 2 q 2 p 2 (k — p) 2 {jp + q) 2 q 2 {p + q) 2 {k + q) 2 

1 



+ 



(p + rf 



+ 



(p + q + r) 



p 2 q 2 p 2 (p + q) 2 q 2 {p + q) 2 J (k 2 + a 2 ) 
For the box diagram we first define 

B(p,q,r) = Bi(p, -q,p + r) . 

Then writing 

B l = B i f ) + a 2 Bf ) + ... 

we have for the scaling part 

Bf\p,q,r)=b\p,q,r) + ^{- I L_ 

1 1 



+ 



I \ 

(r — v) 2 (v — q) 2 J 



r 2 (r — p) 2 (q — r) 2 q 2 (p — q) 2 (q — r) 2 p 2 (r — p) 2 (p — q)" 2 
where b' is the continuum integral 



b'(p,q,r) = / { 

J oo 



k 2 (k — p) 2 {k — q) 2 (k — r) 2 p2g2 r 2 



1 



p 2 (p — q) 2 {p — r) 2 



1 



(k — p) 2 k 2 + a 2 



k 2 k 2 + a 2 J 
+ 2 perms j j 



For the leading scale-breaking piece one gets 

-{ 



B { f\ P , q,r) = ^\ «iB[ 3) (p, <?, r) + (k 2 + ^i) ^ (p, 9, r) ] , 
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with 



Bf{p,q,r) 



where 



and 



£ w (p, q,r) + # w (p,p-q,p- r) 

+ B {i) (q,q-r,q-p)+ B {i) (r, r -p,r - q) , (B.84) 



1 



(p — q) 2 (q — r) 2 (r — p) 2 

{{p ~ q)(p ~ r ) [L(a(q — r)) + c\] + 2 perms} , (B.85) 



cx \ [R(p) - |] 



] = ^\ (p-qn P -r) 2 +2peTmS 



+ 2vr 



Rip) 



IS 



R(k) - | 



(k — p) 2 (k — q) 2 (k — r) 2 



(p — q) 2 (p — r) 2 



1 



(k — p) 2 k 2 + a 2 



+ 2 perms^ } . (B.86) 



Appendix C: 1-loop 4-point functions with O3 or Oq insertions. 

There are many contributions to the 1-loop 4-point functions with insertion of the 
operators O3 or Oq. Some simple contributions are of the form of (vertex or prop- 
agator) correction to the corresponding tree level 4-point functions with operator 
insertion. An other category of simple contributions is where the operators are in- 
serted on external lines of the ordinary 1-loop 4-point function of the model. (Some 
of the contributions belong to both of the above type.) There remain 11 "irre- 
ducible" contributions, which do not belong to either of the above sets of simple 
graphs. They are shown in Fig. [Til Each of these contributions is a complicated 
function of the four momenta and since they are topologically distinct, it is difficult 
to automate the calculation. 

We note that it is sufficient to compute the 4-point functions with insertions 
for special momenta p\ + q\ = = pi + qi- Since some diagrams are singular for 
this configuration one computes with m/0 and then takes the limit m — > (which 
is non-singular) in the final result. We checked that the limits m — > and the limit 
to the special momenta configurations commute, a fact that was (and still is) to us 
not obviously the case. 
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Figure 13: 2-loop contributions to the renormalization of operators 0\ and 02- Full 
circles represent the (4, 6 or 8-leg part) of the operators, full squares stand for the 
(4-point or 6-point) interaction vertex of the model. 
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1 
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0.47746483 


-0.07957747 


0.06287345 


p(C) 





0.09153552 


0.02817918 





^28 


-0.07421002 


0.02872250 


-0.18868881 


-0.07852019 


r 28 


-0.05734358 


-0.10135050 


-0.00290322 


-0.03887267 


^29 


-0.07889581 


-0.11486836 


-0.05199331 


-0.06922970 


^30 


-0.05431565 


0.00882813 


-0.13895289 


-0.07852019 


r 30 


-0.07821415 


-0.07058178 


-0.13051244 


-0.04738831 


^30 





-0.01798155 


0.05047165 






Table 7: Constants for various actions. 
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Figure 14: The 11 irreducible contributions to the 1-loop 4-point function with 
insertion of the operator O3 or Oq. Full circles represent the (2, 4 or 6-leg part) 
of the operators, full squares stand for the 4-point interaction vertex of the model. 
The crossed vertices represent the insertion of 7r 2 , and the starred vertex of diagram 
a represents the insertion of (tt 2 ) 2 . 
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Appendix D: solution of the RG equations (I3.45j) 



In this Appendix we will analyse the solution of the RG equations (|3.45p in detail. 
It is well-known that if there are also integer numbers in the set of differences 
{Aj — A,-} (i ^ j) then in general terms proportional to powers of the logarithm of 
the renormalized coupling may occur in the solution. We will show that nevertheless 
all coefficients we need in our analysis here are actually free of these terms. 

Our starting point is the expansion (|3.46p . where the coefficients k$ may log- 
arithmically depend on the coupling through 

u = In {2p g 2 ) . (D.l) 

In fact, the coefficients are finite polynomials in uj: 

^H=x>r } - CT ' ( d - 2 ) 

cr=0 

where the order k is not larger than the largest integer in the set {Aj — A,}J . Jj. 
Substituting the ansatz (|3.46p into (|3.45p gives the following set of algebraic equa- 
tions for the numerical coefficients: 



(Aj -A t -£ + l)k% a) - (a + l)*f 



1 x CO i (i+l-m) , (m,<x) ( D-3 ) 

m=2 

for £ = 2, 3, . . . and a = 0, 1, . . . , k. (&^' k+1 ^ = by convention.) 

These equations can be solved order by order in the loop expansion. More 
precisely, if we already found the solution up to £ — 1 loop order then in the case of 
Aj — Aj ^ £ — 1 we can find a unique solution for all the ^-loop coefficients kf-^ . If 

Aj — Aj = £ — 1 then we can still solve for all k\^ but k^\ which is not occurring 
in the equations in this case. It remains arbitrary: we will make the solution unique 
by putting it equal to zero in this case. 
We get for the two-loop coefficients 

A, - Aj + 1 : m _ 1 (2 ) ( D - 4 ) 



' 2; A, A, -\) P 



6 Note that this maximal power in lo first occurs at higher loop orders. The actual power of lo is 
always smaller than £, the loop order. 
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and 



A,-A i = l: ,(2,D_ 1(2) (D - 5) 

% ~ 2 ij ' 

The three-loop coefficients are 

^' 2) = o, 



2 : ,(3,D = 1 (2) ,(2,1) ( D - 6 ) 

« 2(A i -A i -2) Pis fl J ' 



1.(3,0) 



1 J fc (3,l) , I J3) , 1 (2) . (2,0) \ (D ?) 

Ai - 2) \*« + 2^' + 2 Pis °i J [ V 



and 



{Aj -Ai-2 
^ 0) =0, 



A i -A i = 2: ~ 2 PiS 2 Pu sj ' ( 8 ^ 

J3,2)_ 1(2)^.(2,1) 
K ij ~ iPis K sj ■ 

Now the general coefficients in ()3.48p depend on oj. This dependence is 
inherited from the w-dependence of the coefficients in (I3.53h . However, the 
coefficients we actually need (the 2-loop coefficients for general n and also the 3- 
loop coefficients for n = 3) are actually free of this dependence. We can see this 
by inspecting the above formulae and remembering the triangular structure of the 
anomalous dimension matrix (and its expansion coefficients). The needed terms are 

iM, \ _ ,(2,0) _ 1 J?) lj\ a) 

which is different from zero for s = 1,2 only and for the case of n = 3 (where we 
need it) 



different from zero only for s = 1, 2 



Appendix E: The Id improved estimator 

It is well known that improved estimators can decrease the statistical error signifi- 
cantly. Hasenbusch [24] proposed an improved estimator for the strip geometry with 
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free boundary conditions in the time direction, which reduces the statistical error 
considerably for the case when £(L) > L, i.e. u = m(L)L < 1. It uses the fact 
that the situation is nearly one-dimensional, i.e. the spins within a time-slice are 
strongly correlated. 

We give here a short overview of the method and add a useful modification used 
in these simulations. 

Consider a new configuration S' where the old spins S(x, t') for t' > t are rotated 
globally, by the same O(n) matrix X(t), while the spins on the time slices t' < t 
are unchanged. Performing such rotations for all time slices, 1 < t < Nt — 1, the 
resulting spin configuration becomes 

S'(x, t) = X(l)X(2) ■■■X(t- l)S{x, t) . (E.l) 

The new spin configuration S' can be described (equivalently but redundantly) by 
the original spins S and the rotation matrices X. 

Consider the MC updates of only the X variables, with fixed S variables. For 
the standard O(n) spin action the action governing the dynamics of X(t) degrees of 
freedom is 

Nt-l 

A X (X\S) =-J2 tr (X(t)Q T (t)) , (E.2) 
t=i 

where the n x n matrix Q(t) is given by the spins of the corresponding two neigh- 
boring time slices, 

Q t] (t) = f3^S l (x,t)S J (x,t + l). (E.3) 

X 

(The interaction terms within a given time slice are O(n) invariant, hence do not 
depend on X and are not written out in (|E.2p .) The improved estimator introduced 
by Hasenbusch relies on the fact that the variables X(t) are independent. Obviously, 

(S'(x, t )S'(y, h)) x = (S(x, t )X(t ) ■ ■ ■ X(h - l)S(y, h)) x 

= S(x, t )(X(t )) x ■ ■ ■ (X(h - l)) x S(y, h) , (E.4) 

where (. . .)x denotes the average over the X(t) variables, each with its own Boltz- 
mann factor 

P(X(t)) cx exp(tr(X(t)Q T (t))) . (E.5) 

Note that updating only the X variables is not ergodic0 - the relative orienta- 
tions of the spins within a given time slice remain unchanged - hence the procedure 
has to be supplemented by other updates, e.g. the cluster algorithm. One does not 
even need to actually update the spin configuration in this way - the main advantage 

7 with the exception of the Id case 
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of the procedure is to yield a very good improved estimator: for Lm(L) < 1 it sup- 
presses the statistical errors considerably better than the cluster improved estimator 

The improved estimator (|E.4h boils down to calculate (or measure) the expec- 
tation values {X(t)) with the weight (|E.5|) . 

(X} = ^ J dXe tviXQT) X, (E.6) 

where dX is the Haar measure over SO(n). It is useful to make a singular value 
decomposition 

Q = UqV , (E.7) 

where U, V £ SO(n) and q is a diagonal matrb|§. Introducing the new integration 
variable Y by X = UYV one gets 

(X) = U(Y)V, (E.8) 

where 



(Y) = J dYe tT ™Y, (E.£ 



For n > 3 this integral is not known analytically, and Hasenbusch [23] proposed to 
measure it stochastically, by an unbiased estimator. 

At this point we slightly modify the procedure proposed in [23]. By making a 
MC update on X, one can start with the value X(t) = 1, since this corresponds 
to the original spin configuration S' = S which was assumed to be an equilibrium 
configuration. Making any number of updates (not necessarily many) starting from 
X(t) = 1 one gets a new configuration S' with the proper weight. The starting value 
Xq = 1 corresponds to 

Y = U T V T . (E.10) 

On the other hand, tr(Yq) depends only on the diagonal elements of Y , hence the 
off diagonal elements in (Y) can be replaced by zero keeping the diagonal elements 
unchanged!! Therefore one can make the replacement 

(X)^U(Y) d .V, (E.ll) 



8 Standard SVD routines return a non-negative diagonal matrix and two O(n) matrices. One has 
to make sure that in this decomposition det U = det V = 1 by changing properly the signs in q if 
needed. 

9 by averaging over properly chosen matrices Y obtained by diagonal SO(n) transformations of 
the form [—1, —1, !,...,!] 
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where Y is an average over the given number of actually performed updates and 
(•••)diag denotes the procedure of replacing the off-diagonal elements by zero. In 
particular, it is possible to avoid any actual updates, just use the estimator 



It turns out that this trick alone reduces considerably the largest eigenvalue of (X) 
and thus the overall stochastic error. 

The Id improved estimator works better as m{L)L decreases, i.e. as the system 
becomes effectively more 1-dimensional. We have compared the errors for different 
estimators in some of our runs at uq ~ 1: 

1: Id improved estimator, exact integration, 

2: Id improved estimator using eq. (|E.12|) with no updates, 

3: Id improved estimator with 100 updates for (X), 

4: cluster improved estimator, 

5: standard estimator. 
The ratios of these errors were: 1 : 1.1 : 2.7 : 4 : 11. This shows that the simple 
procedure described above is quite efficient. 
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